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1 Abstract 

lO ■ In this paper we investigate the limit behavior of the solution to quasi- 

static Biot's equations in thin poroelastic plates as the thickness tends to 
zero. We choose Terzaghi's time corresponding to the plate thickness and 
obtain the strong convergence of the three-dimensional solid displacement, 
fluid pressure and total poroelastic stress to the solution of the new class of 
plate equations. In the new equations the in-plane stretching is described 
by the 2D Navier's linear elasticity equations, with elastic moduli depend- 
ing on Gassmann's and Biot's coefficients. The bending equation is coupled 
with the pressure equation and it contains the bending moment due to the 
variation in pore pressure across the plate thickness. The pressure equation 
is parabolic only in the vertical direction. As additional terms it contains 
the time derivative of the in-plane Laplacean of the vertical deflection of 
the plate and of the the elastic in-plane compression term. 
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1 Introduction 



A plate is a 3D body bounded by two surfaces of small curvature and placed at 
small distance. A plate is said to be thin if the distance between these surfaces, 
called the thickness, is much smaller than a characteristic size of the surrounding 
surfaces. 

The construction of a linear theory for the extensional and flexural deforma- 
tion of plates, starting from the 3D Navier equations of linear elasticity, goes 
back to the 19th century and Kirchhoff 's work. Following a short period of con- 
troversy, the complete theory, nowadays known as the Kirchhoff-Love equations 
for bending of thin elastic plates, was derived. Derivation was undertaken under 
assumptions, which are referred to as Kirchhoff's hypothesis. It reads as follows: 



Kirchhoff's hypothesis: Every straight line in the plate that was originally 
perpendicular to the plate midsurface, remains straight after the strain and per- 
pendicular to the deflected midsurface. 

Since then a theory providing appropriate 2D equations applicable to shell-like 
bodies was developed. Later, due to the considerable difficulties with the deriva- 
tion of plate equations from 3D linear elasticity equations, a direct approach in 
the sense of Truesdell's school of continuum mechanics, using a Cosserat surface, 
was proposed. We refer to the review paper of Naghdi [17] . In classical engineer- 
ing textbooks one finds a formal derivation, based on the Kirchhoff hypothesis. 
For details we refer to Fung's textbook [14]. 

A different approach is to consider the plate equations as an approximation 
to the 3D elasticity equations in a plate domaine fl e = u x (—e,e), where u 
is the middle surface and e is the ration between the plate thickness and its 
longitudinal dimension. Here, for simplicity we suppose that plate is flat and of 
uniform thickness. 

The comparison between 3D model and 2D equations was first performed 
formally in [13]. Then Ciarlet and collaborators developed systematically the 
approach where the vertical variable £3 G (—e,e) was scaled by setting y 3 = 
Xz/e. This change of variables transforms the PDE to a singular perturbation 
problem on a fixed domain. With such approach Ciarlet and Destuynder have 
established the rigorous error estimate between the 3D solution and the Kirchhoff- 
Love solution to the 2D plate equations, in the limit as s — > 0. For details, we 
refer to the article [5], to the book [6] and to the subsequent work for details. 
The complete asymptotic expansion is due to Dauge and Gruais (see [9], [10] 
and subsequent work by Dauge and collaborators). Handling general boundary 
conditions required the boundary layer analysis. For a review, discussing also the 
results by Russian school, we refer to [11]. Further generalizations to nonlinear 
plates and shells exist and were obtained using T— convergence. We do not discuss 
it here and don't undertake to give complete references to the plate theory. 

Many living tissues are fluid-saturated thin bodies like bones, bladders, arter- 
ies and diaphragms and they are interpreted as poroelastic plates or shells. For 
a review of modeling of bones as poroelastic plates we refer to [8]. Furthermore, 
industrial filters are an example of poroelastic plates and shells. 

Our goal is to extend the above mentioned theory to the poroelastic plates. 
These are the plates consisting of a porous material saturated by a viscous fluid. 
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The mathematical model of such a plate consists of the Biot poroelastic equations, 
instead of Navier's elasticity equations. In addition to the phase displacements, 
description of a poroelastic medium requires the pressure field and, consequently, 
an additional PDE. 

Until recently, few articles have addressed plate theory for poroelastic me- 
dia. Nevertheless, in the early paper [4], Biot examined the buckling of a fluid- 
saturated porous slab under axial compression. This can be considered as the first 
study of a poroelastic plate. The goal was to have a model for the buckling of 
porous media, which is simpler than general poroelasticity equations. The model 
was obtained in the context of the thermodynamics of irreversible processes. 

More recently, in [24] Theodorakopoulos and Beskos used Biot's poroelastic 
theory and the Kirchhoff theory assuming thin plates and neglected any in-plane 
motion to obtain purely bending vibrations. A systematic approach to the linear 
poroelastic plate and shell theory was undertaken by Taber in [22] and [23], using 
Kirchhoff's approach (see e. g. [21]). 

In this paper we follow the approach of Ciarlet and Destuynder, as presented 
in the textbook [20] and rigorously develop equations for a poroelastic plate. 

2 Setting of the problem 

We study the deformation and the flow in a poroelastic plate Q = X2, £3) £ 
ul x (— £/2,£/2)}, where the mid-surface ujl is a bounded domain in M 2 with a 
smooth boundary dui 6 G x . For simplicity, we suppose that the poroelastic plate 
Q £ is an isotropic material. E £ (respectively ~E~ e ) is the upper face (respectively 
lower face) of the plate Q e . T e is the lateral boundary, T e = 8ujl x (—£/2, £/2). We 
recall that the ratio between the plate thickness and the characteristic horizontal 
length is e = £/{2L) « 1. 

A poroelastic plate consists of an elastic skeleton (the solid phase) and pores 
saturated by a viscous fluid (the fluid phase). At the pore scale, one deals with a 
complicated fluid-structure problem and in applications we model it using the ef- 
fective medium approach. In poroelasticity the effective modeling goes back to the 
fundamental work by Biot (see [2], [3] and [25]). The deformable porous media, 
saturated by a fluid, are modeled using Biot's diphasic equations for the effective 
solid displacement and the effective pressure. Biot's equations are valid at every 
point of the plate and the averaged phases coexist at every point. The two-scale 
poroelasticity equations were obtained using the two-scale expansions applied to 
the pore fluid-structure equations by Burridge, Keller, Sanchez-Palencia, Auri- 
ault and many others. We refer to the book [19], the review [1] and the references 
therein. Mathematically rigorous justification of the two-scale equations is due 
to Nguetseng [18] and to the papers by Mikelic et al [15], [7] and [12], where it 
was also shown that the two-scale equations are equivalent to the Biot system. 

We note that Biot's diphasic equations describe behavior of the system at so 
called Terzaghi's time scale T = riL 2 c j{kG\ where L c is the characteristic domain 
size, r) is dynamic viscosity, k is permeability and G is the bulk modulus. For the 
list of all parameters see the Table 1. 

In general, for the plate we have two possible choices of time scale: 
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D I NldUL 


DTT A ATTTTV 


Lr 


shear modulus 


V 


drained Poisson ratio 


IG 


inverse of Biot's modulus 


OL 


effective stress coefficient 


K 


permeability 


V 


viscosity 


Li AAlKX £ 


midsurface length and plate width, respectively 


cr — D 1(0 T \ 


small parameter 


T = i]L 2 /(kG) 


characteristic Terzaghi's time 


d 


characteristic displacement 


P = dG/L 


characteristic fluid pressure 


u = (wi,w 2 ,ii 3 ) 


solid phase displacement 


p 


pressure 



Table 1: Parameter and unknowns description 



1. L c = L, leading to T = r/L 2 / (kG) and 

2. L c = £, leading to the Taber-Terzaghi transversal time T ta f, = rj£ 2 /(4:kG). 

If Terzaghi's time is short, then it is necessary to study the vibrations of the 
poroelastic plate. For filters, Terzaghi's time scale doesn't correspond to short 
times, the rescaled acceleration factors max{pj, p s }k 2 G/(r] 2 L c ) are small and the 
acceleration is negligible. Only the time change of the variation of fluid volume 
per unit reference volume ( = ^gP + at div u is not small. For more discussion of 
the scaling we refer to [16]. 

Consequently, we study the simplest model of real applied importance: the 
quasi-static Biot system. 

Following Biot's classical work [2], the governing equations as written in [8] 
take the following form: 

a = 2Ge(u) + div u - ap)I in (1) 

G 

-GAu v divu + «VP = in Q e , (2) 

d k 

-rrilGP + ct div u) Ap = in Q e . (3) 

at 7] 

We impose a given contact force an = V ±e and a given normal flux — — — = U 
at x 3 = ±£/2. 

We recall that the effective displacement of the solid phase is denoted by 
u = (ui,U2,u 3 ), the strain tensor e is given by e(u) = sym y u, a is the stress 
tensor and the effective pore pressure is p. 

Following engineering textbooks approach to Kir chhoff- Love's plate model 
and with some appropriate modifications, we that are able to derive formally 
the poroelastic plate equations. 

For this, in addition to the Kirchhoff hypothesis and following Taber's papers 
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[22] and [23], we suppose that 



the fluid velocity derivatives in the longitudinal direction 

are small compared to the transverse one. (4) 

Following Fung's textbook [14], we perform classical Kirchhoff type formal cal- 
culations and obtain the following equations 

G£A XUX2 u" + g ^ ^ Vs^div^u" + \Zy V ^ N 

2 

+ £(73 + 77V = 0, (5) 

3=1 

a 2 (l-2zA r a(l-2v)£ . u „. 

(7o + 2 G{l-u) )N = \- v div^KX), (6) 

a 2 (l -2v).d , eff N, k d 2 , ff N, \-2vd s 



Gi 3 A2 1 -2i/ a /^ /2 eff , 
^.x^ + A x!,x 2 / /J dx 3 



6(1 -v) X1 ' X2 
„ 2 



e/2 



|E^ + ^)+^ + V, (8) 

where w(xi,X2,t) is the effective transverse displacement of the surface, u w = 

dw 

(ui,Uo), u^(xi,X2,t) — X3— — , j = 1,2, are the effective in-plane solid displace- 

J OXj 

M2 

ments, p e ^ is the effective fluid pressure and N = — p e ^ dx% is the effective 

J-l/2 

stress resultant due to the variation in pore pressure across the plate thickness. 

G£ 3 

We note that D = — is the flexural rigidity of the plate solid skeleton. 

6(1 — v) 

We refer to Appendix for the detailed formal calculation. 

The disadvantage of such approach to derive the equations is in using some 
ad hoc hypothesis. The assumptions (104) that the stresses a^, i = 1,2,3 are 
negligible, can only be satisfied approximatively. They are not even consistent 
with the fact that we use then averages which are non zero. To this classical 
difficulty of the plate theory, we add a new ad hoc assumption (4) on the pressure 
field. 

As in the case of the justification of the elastic plate equations by Ciarlet et al, 
here also the two-scale asymptotic expansion approach gives the correct answer. 
Note that the presence of the effective stress resultant due to the variation in 
pore pressure across the plate thickness, causes stretching even if the boundary 
conditions for m£) are homogeneous. 

Our goal is to extend the Kirchhoff-Love plate justification by Ciarlet et al and 
by Dauge et al to the poroelastic case. Due to its structure, the quasi-static Biot 
system cannot be written as a minimization problem. Therefore, we cannot apply 
the approach of T-convergence to obtain the effective equations. Instead, we start 
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with a weak convergence result, and then, after adding additional correction terms 
show the strong convergence. The obtained result is then of the same type as in 
T-convergence approach i.e. they would coincide in the case of linear elasticity. In 
poroelastic case, we have to deal with more involved norms related to the energy 
of the Biot system. 

In subsection 3.1 we present the dimensionless form of the problem. Subsec- 
tion 3.2 contains a sketch of the proof of existence a unique smooth solution for 
the starting problem. In subsection 3.3 we formulate our convergence results. 
Section 4 is consecrated to the introduction of the rescaled problem, posed on 
the domain Q = u x (—1, 1). Then in Section 5 we study convergence of the 
solutions to the rescaled problem, as e — > 0. In short Sections 6 and 7 we prove 
the strong convergence for the corrected displacement and pressures and stresses, 
respectively. In Appendix 8 we give an ad hoc derivation of the model, which 
follows mechanical engineering textbooks and which is justified a posteriori by 
our rigorous results. 



3 Main Results 

3.1 Dimensionless equations 

We introduce the dimensionless unknowns and variable by setting 
7 = 7 G G; P = —r\ A = — ; T - 



L ' 1-2*/ 4kG' 

Uu £ = u; Pp £ = p; xL = x; tT = t. 

After dropping wiggles the system (l)-(3) becomes 

- A u e - - — v div u e + a V f = in Q £ x (0, T), (9) 

a £ = 2e(u e ) + ( V div u e - ap £ )I in Q £ x (0,T), (10) 
d 

— (7/ + a div u e ) - e 2 Ap e = in Q £ x (0,T), (11) 

where u £ = {u\,u\,u ! ^) denotes the dimensionless displacement field and p £ the 
dimensionless pressure. We study a plate fl £ with thickness 2e = £/L and section 
uj = ujl/L. It is described by 

Vt £ = {(xi,X 2 ,X 3 ) G U X (-£,£)}, 

S +e (respectively S~ e ) is the upper face (respectively the lower face) of the plate 
Q £ . T £ is the lateral boundary, T £ = du x (— e, e). 

We suppose that a given dimensionless traction force is applied on S +£ U £~ £ 
and impose the frictional boundary conditions on T £ : 

a £ n £ = (2e(u e ) - ap £ I + - 2 " ( div u e )J)(±e 3 ) = 

V £ = e 2 (V 1 ,V 2 ,0) + e 3 V 3 e 3 on S +£ ur £ , (12) 
u £ = 0, on r £ . (13) 
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For the pressure p £ , at the lateral boundary T £ = du x (—£,£) we impose a 
given inflow/outflow flux V: 

- \/p £ ■ n = V(x u x 2 ,t) ; n = (711,713, 0) (14) 

and at E +e U S _e , we set 

-|^ = t/ 1 (x 1 ,x 2 ,t) on £ ±£ . (15) 

Finally, we need an initial condition for p £ at t = 0, 

p £ (x 1 ,x 2 ,x 3 ,0) = p in (x 1 ,x 2 ) in ST. (16) 
We write problem (9)-(16) in the variational form: 

LetV £ = {zG if 1 ^ 6 ) 3 ! z| r * =0}.Findu £ G if^O, T, V e ), p e G H^O^T; H\n £ )) 
such that it holds 

/ 2e(u e ) : e{(p)dx + — A div u e div ipdx — a f p £ div <pdx 

Jn s 1 — 2u Jqe Jqs 

= / V e <pds, for every <p G and £ G (0, T), (17) 

Js+ e us- e 



7 / d t p £ (dx + a div d t u £ (dx + e 2 I Vp £ V(dx 



! / K(ds, for every ( 6 H\tt £ ), 
Jan* 



where 



M £ 



V, on r e , 
C/ 1 , on £ ±e . 



p £ \t=o=Pin, in (19) 

3.2 Existence and uniqueness for the ^-problem 

In this section we prove existence and uniqueness of a solution {u £ ,p 6 } G 
H\Q,T\V e ) x H\Q,T',H\&)), p £ \ t=0 = p in of the problem (17) -(19). We 
construct a Galerkin approximation u £ ps v\, p £ ~ Pn and show its convergence 
to {u e ,j9 e } when N — > +00. 

Assumptions 1. p in G # 2 (fi e ), P £ G i? 1 (0, T; iJ 1 (9S7 e \r e )), A/" e G ^(0, T; L 2 (9fi e )). 

Proposition 1. Lei us suppose assumptions 1 . Then problem (17)- (19) has a 
unique solution {u £ ,p £ } G H\0,T;V £ ) x if 1 (0, T; 
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Proof. Let /?,■ be an orthogonal basis of V s with respect to the scalar product 



2e(f) : e(g)dx 



9J 



and £j an orthogonal basis of H 1 ^ 6 ), orthonormal for L 2 (Q e ). Then, the Galerkin 
approximation is 

J 2c(vn) : e{ip)dx + — j div (v\) div ipdx — 

a Pn div (pdx = / V £ cpds, for every ip e V)v, and t G (0, T), (20) 



7 / dtPnCdx + / a div <9 4 (v N )£dx + e / WpN^Cdx 

2 / A/;C^, for every ( G span{£i, ...,6v}, 
Jan* 



p N \t=o =P 



N 



(21) 
(22) 



where pf n is a projection of p in to span{£i, ...,£jv} an d K/v = span{(3i, ...,(3n}. 

Therefore, v N = S^ 1 6 w (t) / 9 z (x) and = S^ 1 s;Ar(t)^(x). Using Korn's 
inequality we obtain that the bilinear form 



A(f , g) 



2e(f) : e(g)dx 



2;/ 



1 - 2u 



div f div sdx 



is ^-elliptic. Let 

A kj = [ (2e(/3 i ) : e(/3 fc ) + div & div (3 k ) dx, l<k,j<3, 

M kj = / £j div p k dx, JC kj =e 2 tte, 1 < j < 3. 

Then, using the definition of and p^, we obtain 



Ab(t) -aMs(t) = ^(t), 

7 -^s(t) + aA< T ^b(t) + Ks(t) = G N , s(0) = s 07V , 
at at 



(23) 
(24) 



where b(t) 



bm{t) 
b NN (t) 



and s(t) 



Siiv(t) 
SNN(t) 



After inserting (23) into (24), we obtain a linear system of ordinary differential 
equations for s = (sijv(t), SNN(t)) T , 

( 7 J + a 2 Af r A~ 1 .A/f)-^;s(t) + /Cs(t) = Q N - aM T A~ l ^ N \t), s(0) = S(W . 

Since the (.A/TA" 1 ./^, £) > c||A4^|| 2 , for all £ e R N , it has a unique solutions 
for all times. Moreover, s e H 2 (0,T) implying the same time regularity for vn 
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and p^. Next, we show a priori estimates, uniform in N. Testing (20) by <9 t v/v, 
and (21) by p^, integrating in space and summing up, we obtain 

9t\ / \e{y N )\ 2 dx + — / | div v N \ 2 dx+ ^ / |pAr| 2 cfel 
Uqc l — lv Jq? z j Qe ) 

+ [ e 2 \Vp N {t)\ 2 dx = d t { [ V £ v N ds-e 2 [ M £ p N ds} 
- / d t V E v N ds+ / e 2 d t N £ p N ds. 

It yields 

|| v Ar||z°°(0,T;V £ ) + || ^ Pn\\l 2 {0,T;L 2 (^)) °°(0,T;L 2 (f^)) 

< C(e) {||P e ||i?i(o,T;i2(an^\re)) + W^Wmio^-^idn^)} ■ 

Next we use (20)-(21) to calculate d^N and d t pN at t = 0. Taking the time 
derivative provides if ^regularity in time. Finally, passing to the limit iV — > oo 
and uniqueness are obvious. 

□ 

3.3 Convergence results 

In the remainder of the paper we make the following assumptions 

Assumptions 2. For simplicity, we assume that p in = 0, that V has a compact 
support in du x (0,T] and that U 1 and V e have compact support in u x (0, T]. 

We start by formulating the scaled limit displacements and pressures: 
Let the scaled in-plane displacements and the pressure {w°, n m } e i^ 1 (0, T; Hq(u)) 3 
be given by 

(7+ 2 1-i/) )7rm+ l-„ div ^KX) = 0, inwx(0,T); (25) 
a / n n\ 1 + ^^ i- / o ox a (l ~~ 2z/) 

= ^(7>i| w =i + ^i|, 3 =-i, P 2 | w =i + P 2 | V8 =_ 1 ), in wx(0,T). (26) 

Theory of the stationary 2D Navier system of linear elasticity yields 

Lemma 1. Under assumptions (l)-(2), problem (25) -(26) has a unique solution 
{w°,7r m }eif 1 (0,T; J ff 1 (u;)) 3 . 

The pressure fluctuation n w G ii 1 (0,T; H 1 ^)) is given by 

(7 + « 2 ^rz^y)^ ~^yj~ a T^7 mA ^ dtW " = 0inQx (°' T )' ( 2? ) 

9y 3 ^w\ y3 =i = dy 3 7T w \ y3= -i = -U 1 in (0,T), 7rJ i=0 = in Q. (28) 
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and the vertical displacement w® e if 1 (0,T; Hq{uj)) by 



/ S1 , K W 3 + 1 _ ~ A 2/i,2/ 2 y i yavr dy 3 = P 3 | y3 =i + ^3^3=1+ 
diyy 1 ,y 2 (n\y 3 =i + V 1 \ y3= - 1 ,V 2 \ y3=1 + V 2 \ ys=1 ) in Ox(0,T), (29) 

Proposition 2. Under assumptions (l)-(2), problem (27)- (29) has a unique 
solution {w% tt w } G H^O, T; # 2 (u;)) x tf^O, T; Lg(O)), 5^7r M G ^(0, T; L 2 (fi)) ; 
/_ x 7r^ dy 3 = 0. 

Proof. See Proposition 10. □ 

Theorem 1. (Convergence of the displacement and the pressure,) Let 

{w°,7r m } be given by (25) -(26) and {w^,ir w } by (21) -(29). Furthermore let 

no / \ a(l — 2u) X3 . , 

E cor (Xi, X 2 , X 3 , t) = —r r (7T^(Xi, 2 2 , — , t) + 7T m (Xl, X 2 , t))- 

2(1 — v) e 
y3^( div xl;X2 w°(x 1 , x 2 , t) - —A XuX2 w%(xi, x 2 , t)). (30) 

Then we have the following convergences 

1 f du £ 

lim max — — / e~ 2 \ — ± - eE°\ 2 dx = 0, (31) 

E -H)0<t<T I fie J J Qe 1 Qx z 

lim max — — / ( — - w° + — V XUX2 w° J | 2 = 0, 1 < i, j < 2, (32) 

1 f | e j3(u £ ) |2 



lim max — — / I J v ; I 2 dx = 0, 1 < j < 2, (33) 

£ ^oo<t<T |fie| y n£ 1 e 1 -j- v / 

lim max -j- / 1^ - w° + — V Xl X2 w°| 2 dx = 0, 1< j < 2, (34) 

& <t<rjjF| j[ l M 3- ^s^i'^^)! 2 dx = 0, (35) 



lim max — — / I ir w (xi,x 2 , — , t) — n m (xi, x 2 , t)\ dx = 0, (36) 

e^oo<i<r \\i £ \ j ne e e 



lin hB3r TriU / {l V ^i-^P £ | 2 + \ d x 3 P £ - d y:i n w \ y3=X3/£ \ 2 } dx = 0. (37) 

Proof. It is a direct consequence of Proposition 12 . It is enough to rescale back 
the variables and the unknowns. □ 

Remark 1. Convergence (31) justifies Taber's hypothesis (4). 

Theorem 2. ('Convergence of Biot's stress,) Let us suppose assumptions 
from Theorem 1. Then we have 



— / ^ - 2e 3j w° + 2-±— ^ 3 
e ^oo<t<T |fi £ | J n , e e dx) 



lim max ^ / - 2e ij (W > ) + 2——^ + 2E° cor \ 2 dx = 0, j = 1, 2, (38) 



lim max-L/ |^ - 2e 12 (w°) + 2^ 7 ^-| 2 dx = 0, (39) 



e^oo<t<r \il e \ J Q£ e e dx\dx 2 

lim max JL / \°1 



1 f o~ £ 

lim max — — / |^| 2 dx = 0, j = 1, 2, 3. (40) 
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Proof. It is a direct consequence of Proposition 13 . It is enough to rescale back 
the variables and the stresses. □ 



Remark 2. We note that convergences (40) justify Kirchoff's hypothesis for the 
case of the poroelastic plate. 



4 The scaled problem 

Our objective is to study the behavior of the displacement u e and the pressure p e 
in the limit e — > 0. Since the problem is defined on the domain Q £ = 00 x (—£, e), 
which is changing with e, it is convenient to transform it to a problem defined on 
a fixed Q = 00 x (—1, 1). 

To the unknowns u £ and p £ , defined on Q £ , we associate the scaled displace- 
ment and pressure fields w(e) and ir(e) by the scalings yj = Xj for j = 1,2 and 
2/3 = if w ith x G Q £ , y G fl such that 

Wj{€){y 1 ,y 2 ,y3,t) = , for j = 1,2, 

«>3(e)(j/i,3/2,2/3,*) = w|(a:i,X2,X3,t), 

Tr(e)(y 1 ,y2,y3,t) = . (41) 

By direct calculation we find out that the scaled displacement and pressure fields 
satisfy the following variational problem: 

Let V(Q) = {z G H 1 ^) 3 ] z| r = 0}, where T = du x (-1,1) and denote 
w = (wx, w 2 ) and <p = (</?i, ip 2 )- 

Find w(e) G ^(0,7,^(0)), vr(e) G H\0,T; H^Q)), such that it holds 

2 / % ( w ( £ ))%(^)^ + 1 — ^ / div !/i,2/ 2 w(£) div w>M ^dj/ 
iJ=1 1 - ^ Jn 

e) divj, liW ^dy+ / S e 3j (w(e))e 3j (^)dy 

-« / K(£)d y3 (p 3 dy + — / dw yu y 2 w(£)dy 3 (p 3 dy 

Jn v — iv Jq 

+ 1~~2^ Sn dysW3 ^ div y^ dy } + (1 - 2i/)e 4 ^ dy3W3 ^ dy3<Psdy 
2 

V f Vjtpjds + / PaVada Vy? G V,t G (0,T), (42) 

7 / d t Tr(e)(dy+ / a diVj /1 ^ 2 ^w(e)C^ + e 2 j V yuy2 7r(e)V yuy2 (dy 
Jn Jn Jn 

+ae~ 2 J d V3 d t w 3 (e)( dy + J — — dy = -e J J V(ds 

±U 1 ( ds, VC G if 1 (ft), (43) 
7r(e)|t= = 0, in fi. (44) 



-a n 
In 
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5 Convergence of the scaled displacement and 
the pressure as e — > 

Proposition 3. The following a priori estimates hold 



H w ( £ )||if 1 (o. r ;-H' :l (") 3 ) 


< c, 




(45) 


\\ 7r ( E )\\H 1 (0,T;L 2 (n)) 


< c, 




(46) 


\hj(w(e))\\ L 2 {nx{0iT)) 


< c, 


l<i,]<2, 


(47) 


ll e i3(w(e))|| L 2 (Qx(0iT)) 


< Ce, 


1 < i < 2, 


(48) 


||e33(w(e))||i3 ( nx(o,T)) 


< Ce 2 , 




(49) 


\\dy 3 n(e)\\ L 2 {nx{0tT)) 


< c, 




(50) 


ll v w,io 7r ( e )lli- 3 (nx(o,T)) 


< £. 

e 




(51) 


equations (42) and (43) by <p = d t w 


^e) and ( = ir(e) 


, respectively, 



integrating them and summing up, we obtain 

1 d 
2dt 



{ 2 J lb M w (^))| 2 */+ Y~~2^ J I diy yi^( e )\ 2 d y + 7 J \n{e)\ 2 dy 



e 2 (l-2u) 



div yi ,y 2 w(e)dy 3 u 3 (e)dy 



\d y3 w 3 (e)\ 2 dy\+e 2 / |V 



2(1-") 
(1 - ^)e A jn 

tAYI I Vj w j(£)ds + / V 3 w 3 (e) ds 



/ 


9tt(£) 




dy 3 



dy 



-E 



dtVjWj(e)ds — / d t V 3 w 3 (e) ds 

E+US- is+US" 



if/^e) - £ / / yvr(e)rfs. 
s+us- J-i iaw 



We can estimate the terms on the right hand-side of (52): 





(52) 



±U 1 n(e)ds 



s+us- 



n 

/ V 3 w 3 (e)ds = d y3 (w 3 ( 



{U x n(e))ds 



< C 



_d_ 



-7T 



(53) 



I 2/3 + 



dy 



< 



C / (\d y3 w 3 (e)\ + \w 3 (e)\) dy < C( / |e(w(e))| 2 dy) 1 ' 2 , 

Jn Jn 

using Korn s inequality 



(54) 
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d t V jWj (e) ds = | / d y3 ( Vj{y3=1 V ^=~ 1 
s+us- Jn 

^3 1 2/3=1 + "P?'ly3=-i 



f) I ' 3 13/3 = 1 ' j|y 3 = -l 

"ya I 2 ^ 3 



Wi(e)tZj/| < C / \d y3 Wj{e)\ dy + C / \w 3 (e)\dy < 



n 

1/2 

|2 



C U \V Wj {e)\ dy \ < C\\e(w(e))\\ L2{n) . (55) 

The initial value of w(e) is calculated using p in and we have w(e)\ t=0 = 0. Using 
equation (43) we find that d t 7i(e)\ t= o = 0, as well. The integration in time of (52) 
gives estimates (45)-(51), but with L 2 -time norms. 

Next we calculate the time derivative of equations (42)-(43), test by (f = 
d a w(s) and ( = d t TT(e). After repeating the above calculations, we obtain esti- 
mates (45)- (51). □ 

Proposition 4. The weak limit {w*,tt } G H l (tt x (0,T)) 3 x H l (Q, T; L 2 (Vt)), 
dy 3 ir° G i/ 1 (0, T; L 2 (Q)) of the weakly converging subsequence {w(e),n(e)} £> o 
belongs to H 1 ^, T; Vm(£1)) x H\0,T; L 2 (n)), where 

14 1 (fl) = {vGff 1 (fl) 3 ,e j3 (v)=0 in Q and v = on dcu x (-1,1)}. 

Proof. Using weak compactness, we find out that there exists a subsequence, 
denoted by the same subscript e, and elements w* G H 1 ^ x (0,T)) 3 and 7r° G 
H 1 ^ x (0,T)) such that 

w(e) w* weakly in H\Q x (0, T)) 3 as e ->• 0, (56) 

?r(e) 7T° weakly in if^O, T; L 2 (tt)) as £ 0, (57) 

d y3 n(e) -± d y3 n° weakly in H 1 ^, T; L 2 (tt)) as e -> 0, (58) 

e(w(e)) e(w*) weakly in L 2 (fi x (0, T)) as e 0. (59) 

Using (48)- (49) we obtain 

ll e 33( w (^))IU 2 (^x(0,T)) < Ce, J = 1,2, (60) 

||e33(w(e))||i2(nx(o,T)) < C^ 2 - (61) 
Hence, by weak lower semicontinuity of norm 

||e j3 (w*)|| L 2 (nx(0iT)) < limiiif ||e i 3(w(e))|| L2 (nx(o,r)) = 0. (62) 

Consequently, w* e H\0, T; V k i{Q)). □ 

Lemma 2. The space Vki defined in Proposition 4 is characterized by the following 
properties 

V kl (ty = {v G H\Q) 3 , v 3 = g 3 -y 3 ^, 9j G Hfo), j = 1,2; v 3 = g 3 G # 2 (u/)}. 
Proof. The proof of this Lemma is given by Ciarlet [6] [p. 23]. □ 
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Corollary 1. There is w° G H\0,T; H^(u) 3 ), w° G H 1 ^, T; H^(u)), such that 
inQx (0,T) 

w j =Wj(yi,V2,t) -y 3 d y wl(y 1 ,y 2 ,t), j = l,2; w* 3 = wl(y 1 ,y 2 ,t). (63) 



Proposition 5. Let < n° >= ~ f_. it dy 3 . Let 



E CO r = E cor (y 1 ,y 2 ,t) 



a(l - 2v) < 7T° > -2vdiv 



y\,yi^\ w \i w 2) 



2(1-1/) 
Then, 

1 <9w 3 (^) 1 «(1 - 2z/)tt° - 2z/ cfe W)1/a w* 
^e 33 (w(,)) = - 2{i _ u) = E cor = 

E C or + ^A yuy2 w° + ^^(tt - < vr >), (64) 
weakly in L 2 (Q x (0,T)), as e — >• 0. 

Proof. After setting <y2j = in equation (42) and multiplying by e 2 , we get 



2 / ^2 e 3j( w ( e )) d y 3 i P3 dy-a 7r(s)d y3 (p 3 dy 
Jq j =1 Jn 

^ f JS „ ^« ,„ J .. , 2(1-1/ 



^ div w , w w(e)9 OT ¥? 3 + ^ J ^d y3 w 3 (e)d y3 (p 3 dy 



1 - 2f 

e 2 / PaVads. 



js+us- 

Passing to the limit £ — )■ yields 

' 2(1 -v) 2v \^ 3 




o Jn 



1 _ iv E * cor + 1 - 2v diVvi ^ W *' W ^ ~ a7T ) ~dy 3 ~ dy = ° 



for all <f 3 G C(0, T; if 1 (fi)), </? 3 | w = 0. Now it is straightforward to conclude 
(64). □ 

Proposition 6. It holds 

ej3(w(e)) _^ ^ wea Hy in L 2 (Q x (0,T)), as e -»• 0, j = 1, 2. (65) 
Proof. Setting y? 3 = in (42) and multiplying by £ yield 



In the limit £ — )■ we obtain 

eia(w(e)) 



X i; weakly in L 2 (fi x (0,T)) (66) 



and J n Xj9 y3 ^Pj = 0, Vy>j G if 1 (fi), v 9 ileajx(-i,i) = 0. It finishes the proof of 
Proposition 6. □ 
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Proposition 7. (w^w®) satisfies system (25)-(26). 

Proof. We test equation (42) by <p G Vkl(Q)- Since ej 3 (cp) = 0, we have 

2 / e iA w ( £ )) e M d y + i 2 *1 / div yi,j/ 2 w(e) div yii?/2 ^ 
./n iJ=1 1 - ^ ./n 

-a / ttO) div^^ tfy + 2 ^ / d y3 w 3 (e) dw yuy2 fi dy 
Jn K 1 ~ Av ) e Jn 

2 

• =1 Je+us- Je+us- 
We use Proposition 5 and pass to the limit e — > in equation (67). It yields 



i div yiiJ/2 w*) div^^dy = V] / P^j ds+ / P 3 </?3 <is (68) 
1 _ ^ " Je+ue- Je+ue- 

and choice <p = (<?i,g2,0), gj G i/g(w), gives equation (26). 

Finally, we take ( = C(2/i? 2/2, ^) as a test function in (43). Using Proposition 
5 and zero initial data, we obtain the equality (25). □ 

Proposition 8. The pressure equation reads 

or/ 2 1 _ 2zy \ n l-2v , n nN , d 2 7T° 

1 - 2i/ 



2/3^,^,^ = in Ox(0,T), (69) 

^ s 7r | ttJ=1 = ^ s 7r°| lft , = _i = -i7 1 (j/i,j/2,t) m (0,T), (70) 
ttVo = m n. (71) 



Proof. Passing to the limit e — > in equation (43) yields 

7 / d t TT°Cdy+ / a( div^^w* + dtE^Qdy 



+ / = - / ±U 1 (ds, VC G ff^fJ). 

Jn dy 3 oy 3 Je+ue- 



(72) 



Using Proposition 5 and zero initial data, we obtain from (72) system (69)-(71). 

□ 

Corollary 2. The function n w = tt°— < tt° > satisfies system (21) -(28). 
Proposition 9. The limit plate deflection w\ satisfies equation (29). 
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Proof. We take as test <p = (-y3— -, -2/3-^,03), with 5-3 G H 2 (uj). It yields 

/ \ d 2 g 3 d 2 g 3 d 2 g 3 

enlV 5 ) = -y3-aT2» e 22(^J = -2/3-^-2-, ei 2 (v?) = -y 3 



div ^1,2/2(^1^2) = -y 3 A yuV2 g 3 



Passing to the limit e — > in equation (42) yields 



2 

f fo 2XT 92W 3 d2 9 3 , !- 2l/ / . 2z/ ?/3 a On a \ J 



Z )J — 1 

2 

= -J2f V j y 3 ^ 1 ds+ I Psgsds. (73) 
~[ Js+us- t% Je+us- 

Equation (73) implies (29). □ 
Proposition 10. System (21) -(29) has a unique solution. 

Proof. It is enough to study the problem with the homogeneous data. Let 
(u° 3 ,p°) G H 1 {Q,T]H 2 (uj)) x H\Q,T;H\Sl)), < p° >= 0, be a solution. We 
take g 3 = d t u 3 in (73) and use £ = p° as a test function for system (27)-(28). 
After summing up we obtain the equality 

d t [ -^-\A yi , y2 uT dy + [ ^(1 ~ 2 ») p o Ayi yAu o dy+ f \ dmP y dy+ 
Jn 1 ~ v Jn 1 — v Jn 

\ dt i {l + a2 W^V) ){p0)2 dy " L y -^T^ p ° A ^ Au " dy = a (74) 



After noting that the second and the fifth terms cancel each other, we conclude 
that p° = and A yi y2 u\ = in to. The last equation, together with u 3 G Hq(uj), 

4 



yields u9 = 0. □ 



Proposition 11. The whole sequence {w(e),p(e)} satisfies 

Wj(e) -± w° - 2/3—^-, j = 1, 2, weakly in H 1 ^ x (0,T)) as e -»■ 0, (75) 

iw 3 (e) w 3 ° wea£% m if 1 ^ x (0, T)) as 5 -> 0, (76) 

7r(e) 7T° = vr m + tt w wea£% in if^O, T; L 2 (fi)) as e -> 0, (77) 

dya^e) 9 y3 7r° weakly in H\0,T; L 2 {Q)) as e -»■ 0, (78) 

where w° and ir m =< tt° > are given by by (25) -(26) and {w 3 ,ir w } fry (^27)- 



6 Strong convergence 

In this section we establish that the weak convergences from previous section 
imply the strong convergences. 
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With such aim, we introduce the corrected unknowns, for which the weak 
convergence to zero was already established in Subsection 5. 

Let ty £ G Cq°(u;) be a regularized truncation of the indicator function l w , 
equal to l u if dist (y,du) > e and such that 1 1 ^ yi ,y 2 ^ e \ = Ce 1 ^" 1 . We set 

^(e) = Wj (e) - w Q j + y 3 -^p, J = 1, 2; 

Ue) = w 3 (e)-w° 3 -e 2 S, J (79) 
«(e) = 7r(e) — 7r°, 

with £ = ^ £ (yi,y 2 ) / E* or (y 1 ,y 2 , a,t) da. The choice of correcting terms is 

Jo 

explained by the following result 
Lemma 3. We have 

-e j3 (£(») wea£% in fl^O, T; L 2 (fi)), j = 1, 2, 3, as e ^ 0, (80) 

£(V)^0 weakly in H l (0, T; ifg(fi)) 3 , as £ ->• 0, (81) 

^e 33 (^(£))^0 weakly in H\0, T; L 2 (Q)), as e -> 0, (82) 
Proof. First, using definition (79) we find out that 

-ejMe)) = -ej 3 (w(e)) - 2^7^> ' 

J by(65) 

implying (80). 

(81) follows from Proposition 4 and Corollary 1 and smallness of the correcting 
terms. 
Finally, 

^e 33 (^)) = ^33(w(e))-^^0. 



□ 



Next, using equation (68) we obtain 

2 / J^e^w -y 3 Vj, li j, 2 ^)e ii (^)dy + — ^— / div^ li?/2 (w 
./n w=1 l - ^ Jn 

-yzV yij y 2 wl) div yii y 2 $dy - a / tt° div yuy2 $dy + 2 / V — e 3j (cp)dy 

Jn Jn j=1 °Vj 

+ \{~ a I ^dy^sdy + 2U I div yi> y 2 (w° - y 3 V yi) y 2 w$)dy 3 ip 3 dy 
e k Jn 1 ~ Zl/ Jn 



fd(w^ + e 2 S) , _,i 2(l-i/) /'^ . n 2(na 

7n — ^ — dlVyi ' y2 ^ + (i - 2z/)£ 4 / Q 2/3 + 5 £ > d y^ dy 

V / + / Ps^ds + y / (tf e - l)£* or div^,^ 

~[ Je+us- Js+us- 1 — lv Jn 



1 - 2z/ 

2 
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f d£ f 1 2z^ f 

+ 2 1 E dy-^J Q -°^3 rfy + ?T -^ ^ div,^^ - 

ysV^wg)^^ + ( fl^)! 2 J^AAmdy, V<peV,te(0,T). (83) 

Consequently, we observe that {£(e), ^(e)} satisfy the following variational equa- 
tion 

f 2z^ f 

1 /" 2 



—a / «(e) div^^^df/ 
'n 



-a J K(e)dy. A (p 3 dy + - y div^^^e)^^^ 

+ 1~~W X div ^ 2< ^ y } + (i _ 2v)e^ J Q d y&( £ ) 9 y3<P3dy = 

f d£ ot f 2z^ f 

2 / E d^ e ^ dy+ i n ° d y^ dy - — — / (W e - l)£ c * or div^p cfy 

*S & l j \ 3 wis */ 16 



2(1-1/) 



/ VeE^dntpady, V<peV,te (0, T). (84) 



(1 - 2u)e 2 
Effective pressure equation (72) implies 



7 f d t n°C dy+ [ a{ dw yi , y2 d t w* + e- 2 d t d y:i (w° 3 +e 2 £))(dy 
Jn Jn 



q <9*/3 ^2/3 ^ + £ 7n V ^ 2?r V w.»^ dy = e V yit y 2 7r V yu y 2 ( dy 

- [ ±U\ds + a f(* e - l)d t E: or ( dy, V( e H\Q). (85) 
je+us- Jn 

Consequently the pressure equation for ^(e)} is 

7 / d t K(e)(dy+ / a <tiv yi;y2 d t i{e)Cdy + e 2 / V ' yi , y2 n{e)V ' yi , y2 (dy 
Jn Jn Jn 

+ae ~ 2 L d » d <^ )c * + L * = - e L L nds 

-e 2 [ V, lty yV, ltK C dy - a /(*« - l)cMCC dy, V( 6 »'(!)). (86) 

if) if] 

By the choice of the correcting terms, £(e) =0 on du x (— 1, 1). Now we test (84) 
by <f = dt£(e), (86) by £ = and add the obtained equalities, to obtain the 
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following energy equality 

l_d_ - 2 
2dt 



{2/EWfW)P* + T ^/ 



+- 



X>3;(£(e))l 2 dy + 



2(1-^) 
(1 - 2v)e 4 

2v 



(1 - 2z/)e 2 



div 2/i,2/ 2 £( £ ) dy + -f \k(e)\ dy 
Jn 

div yi ,y 2 £{£)dy 3 £ 3 {e) dy 



\dy 3 ^ 3 (£)\ 2 dy\ +e 2 / \V yi , y2 n(e)\ 2 dy 



1 


dK,(e) 


In 


dy 3 



dy 



1 - 2v 

2(1-") 
(1 - 2v)e 1 



div. 



£ Jo, 

(* e - l)£ c * or div yi;y2 f(e) ) dy 



a 
2v 



d vs&( £ ) dy 



2(1-1/) 



1 - 2i/ 



(1 - 2v)e 1 
d t div 7/1 ,„fw u - VzV yi ,y 2 w 



v e d t E: or d y M £ )dy + 



(tf e - l)9^: or div^e) dy 



/ V W)W 7r V Wj , /3 K(£) dy-e / yft(e)ds - a / (* e - l)<9 t £* or .K(» dy. 



Proposition 12. For the whole sequence {C,(e),k(e)} we have 

e ;i(£( e ))^° strongly in tf^O, T; L 2 (n)), i,j = 1, 2, as e ->• 0, 
£(e)^0 strongly in H\0, T; L 2 (Q)) 3 , as E -> 0, 

ea(£(<0) 



e 2 <9y 3 



^0 strongly in H 1 ^, T; L 2 (fi)), as e^O, 
^0 strongly in H 1 ^, T; L 2 (fi)), as e^O, 



(88) 
(89) 

(90) 
(91) 



/c(e)->-0 strongly in H x (0, T; L 2 (fi)), as e^0, (92) 
w «(e)-»O strongly m H\0,T; L 2 {{})), as e ->■ 0, (93) 
eV Wjya K(e) -»• strongly in H 1 ^, T; L 2 (fi)), as £ ->• 0. (94) 
Proof. We follow proof of Proposition 3 and use equality (87). In fact only new 

term to be estimated is e / / ^K(e)ds. We recall the well-known interpolation 



inequality 
It yields 



||ClU 2 (0a;x(-l,l)) < C||C|li2 2 (Q)| 



1/2 




V7c(e) ds| < Ce||^||^ (a-x( _ lil)) ||K(e)||^ (0) (||«(e)||^ n) + 



1 1/2 



i Ja 



11^,^)11^) < C^Wk&Ww) + yl|V WfIft K(e)llzw- (95) 



57) 
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Now we integrate equality (87) in time, use Proposition 11 and Lemma 3 and 
conclude the strong convergence in L 2 (Q x (0, T)). Iterating the argument after 
calculating the time derivative of equations (84) and (86), yields the convergences 
(88)-(94). □ 



7 Convergence of the poroelastic stress 

The rescaled stress er(w(e)) is, in analogy with (10), given by 



a(w(e)) = 2e(w(e)) + 



2v 



1 - 2v 



div w(e) - aix(e))I in Q x (0,T). (96) 



Nevertheless, this quantity does not correspond to the rescaled dimensionless 
physical stress a 6 . We introduce the rescaled poroelastic stress a{e) by 



2^M + D{e) 2e 12 (w( £) ) 

oyi 



i 



ei 3 (w(e)) 



e 

2e la (w(e)) 2^^ + J D(e) -e 23 (w(e)) 

-e 13 w e -e 23 w £ -2"-^ + ^(e 

£ e e 2 dy 3 



(97) 



where 



D(£) 



2v 



div2 /liy2 w(£) - an(e) + 



2u 1 ^W 3 (£) 



(98) 



1 2// 1 ' ' l-2^£ 2 dy 3 ' 

As a direct consequence of Proposition 12, we obtain the following conver- 
gences for the stresses 

Proposition 13. Let {w°,7r m } be given by (25)-(26) and {w%ii w } by (27)- (29). 
Let E* or be defined by (64)- Then we have 



D(e)^-2E* cor in H\0, T; L 2 (n)) as e -> 0, 



— 2e 3 -j(w J 



o 

3 1 2E* cor -+ 



2 ^ 3 a 2 



in H l (Q,T;L 2 (ti)) as e^O, j = 1,2, 



012(e) 



2e 12 (w°) + %3 



- 3 ->0 in i/VO^L 2 ^)) as £ ->■ 0, 



o-i3(e) 



->■ in tf^O, T; L 2 (fi)) as £ -> 0, j = 1, 2, 



033 (£) 



->0 in H\0,T;L 2 (Q)) as e -> 0. 



(99) 

(100) 
(101) 

(102) 
(103) 



8 Appendix: A classical Kirchhoff type plate 
equations derivation 

We follow [14] and, together with Kirchhoff's hypothesis, suppose that 
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1. The vertical deflection of the plate takes form u 3 = w = w(xi,X2,t). 

dw 

2. For small deflections w and rotations , ^2) we have $1 = — — and t? 2 — 

ox 2 

dw 

. Then, Kirchhoff 's hypothesis implies that the displacements u satisfy 



dx\ 



dw dw 
u x = u"(x l7 x 2 , t) - x 3 — ; u 2 = v%(xx,x 2 , t) - x 3 —, (104) 

OX\ OX2 

where (u^u^) are the tangential displacements of points lying on the mid- 
surface. Now a direct calculation gives e^u) = 0, j = 1,2. Since 

a = 2Ge(u) + (- — — div u — ap) I, we conclude that also Uj 3 = 0, j = 1, 2. 

3. For the pressure field, we impose the normal velocities U e at the top and bot- 
tom surfaces and suppose (4). Finally, 033 is also supposed small throughout 
the plate. The later assumption gives 

e 33 (u) = -(a u + 022) + —^77 tP- (105) 

6iK ' 1 + ^2G V ; 2G(1 + z/T 

4. For the other components of the strain tensor we have 

_ d< <9 2 w 1 <9< <9u£ <9 2 «7 

611 U dxi Xs <9x 2 ' 612 U 2 <9x 2 ctei Xs 8x18x2' 

. . <9-u 2 <9 2 u; . . f . . , ... a(l - 2v) 

e 22 (u) = _ - x 3 ^ ; e 3 3(u) = - — (e U (u) + e 22 (u)) + ^-^yP; 

2z/G , 2Gv . . , , a(l - 2i/) 
div u - = (en(u) + e 22 (u)) p 



l — 2i/ 1 - z/ v v ' v " 1 - v 

5. For the other components of the stress tensor we have 

2G , , , . a(l-2v) . . 

a n = - (en u + ^e 22 u - -\ >-p 106 

1 — v 1 — v 

^22 = 7 e 22 u + i/en(u)) - ^ 107 

1 — v 1 — v 

8n u 8n u 8 2 m 
a l2 = 2Ge 12 (u) = G(p- + p.) - 2Gx 3 ^-. (108) 

<7x 2 axi 8x18x2 



Next we define the stress resultants (forces per unit length) N\, iV 2 , AT 12 by 

-£/2 J-^/2 



an cfe 3 , N 2 = a 22 dx 3 , N12 = / cr l2 dx 3 , (109) 

-£/2 7-^/2 J-t/2 



the effective stress resultant due to the variation in pore pressure across 
the plate thickness N by 

rt/2 

N = - pdx 3 (110) 
J-e/2 
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and the external loading tangential to the plate 

fi = 0- i3 \x 3 =l/2 - <?i?,\x 3 =-l/2 = Vi + , 2 = 1,2. (Ill) 

Next we average the equation (2) over the thickness and obtain 

^ + f!i + /l = 0, (112) 
OX % ox 2 

dN 12 dN 2 , . . 

ifef + i^ + A = ' < 113 > 

Inserting (106)-(108) into the formulas (109) gives 

2G£ .duV du%. a(l-2i/) Ar A . 
1 — V OXi ox 2 1 — v 

N2 = + v^-) + A. L N; 115 

\ — V OX 2 OXi 1 — V 

A substitution of equations (114)-(116) into the (112)-(113) yields the equations 
for stretching of a plate of uniform thickness (5). 

We need one more equation to complete the system (5). We have 

1 — 2v , . . a(l — 2v) ..... 

U = 1 _ y 1 dn W 2 K> W 2) -XgA^^Wj + — - _ p, (117) 

div u rfx 3 = * { ] - 2U) div^^K,^) - ^!~ 2 ^ iV. (118) 

1 — 1/ ZLr^l — V) 

Next we average equation (3) over thickness, use the assumption (4) and expres- 
sion (118) and get (6). Inserting (117) in equation (3) and using hypothesis (4) 
and equation (6), yields equation (7). 

It remains to find the equation for transverse deflection and for the bending 
moment due to the variation in pore pressure across the plate thickness. 

The bending moment M due to the variation in pore pressure across the plate 
thickness is given by 

M2 

M — — x 3 p dx 3 . (119) 
J -1/2 

For the stress moments of the plate, which have as physical dimension the moment 
per unit length, we have: 

1. The twisting moment M\ 2 '■ 

f e/2 G£ 3 d 2 w 
Mi 2 = / a 12 x 3 dx 3 = — . (120) 

J -1/2 13 OX\OX2 
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2. The bending moments Mi and M 2 : 

f e/2 , Ge ,d 2 w d 2 w, l-2i/ w , . 

Mi = l /2 ^ ^ = -w^m + * &?> + a — M (i2i) 
= l /2 **** = -6(1310 ( m + ^ + a — M - ^ 

3. The transverse shear {Q\,Q-2): Q% = fl^/ 2 o~i3 dx^, i = 1,2. 

4. Resultant external moment (m 1 ,m 2 ) : m ; = + )■ 

Next we multiply the equation (2) by X3 and integrate over the thickness to obtain 

™L + °*>- Ql + mi = 0, (123) 

OXi OX2 

dM l2 dM 2 „ , N 

+ - Q2 + m 2 = 0. 124 
ctei OT 2 



Averaging the the third equation in (2) over the thickness yields: 

V e 3 + P~ e = (125) 



dQi , <9Q 2 



<9xi <9x 2 

Following Fung's textbook, we eliminate Qi from (123)-(125) and obtain the 
equation of equilibrium in moments: 

d 2 M 1 d 2 M l2 d 2 M 2 dm x dm 2 w m , , , 

ctef ax\ox 2 ox 2 ox\ ox 2 

After inserting the formulas (120)-(122) into (126), we obtain the poroelastic 
plate bending equation (8). 
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